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A   NEAR-FIELD   ARRAY   OF   YAGI-UDA  ANTENNAS 
FOR   ELECTROMAGNETIC   SUSCEPTIBILITY  TESTING 

David  A.   Hill    and  Galen  H.   Koepke 

In  electromagnetic  susceptibility  testing  of  electronic  equipment, 
the  ideal  incident  field  is  a  plane  wave.  To  approximate  this 
condition,  a  seven-element  array  of  Yagi-Uda  antennas  has  been 
constructed  and  tested  at  a  frequency  of  500  MHz.  The  element 
weightings  are  determined  by  a  near-field  synthesis  technique 
which  optimizes  the  uniformity  of  the  field  throughout  a  rectangu- 
lar test  volume  in  the  near  field  of  the  array.  The  amplitude  and 
phase  of  the  electric  field  have  been  measured  throughout  the  test 
volume  with  a  short  dipole  probe,  and  the  agreement  with  the 
theory   is  excellent. 

Key  words:  antenna  array;  array  synthesis;  dipole  antenna; 
electromagnetic  susceptibility  testing;  mutual  coupling;  near 
field;   plane  wave;    Yagi-Uda  antenna 


1.      Introduction 

In  electromagnetic  susceptibility  testing  of  electronic  equipment,  the 
ideal  incident  field  is  a  plane  wave.  Far-field  antenna  ranges  are  generally 
unacceptable  because  they  waste  power  outside  the  test  volume  and  because  the 
far-field  distance  is  often  inconveniently  large.  An  ideal  test  facility 
would  supply  most  of  the  power  to  the  test  volume  containing  the  equipment 
under  test  (EUT),  and  the  test  facility  dimensions  should  not  be  too  much 
larger  than  the  test   volume  so  that   it  could  be   located  indoors. 

The  compact  range  [1]  uses  a  paraboloidal  reflector  to  produce  an  ap- 
proximately plane  wave  in  the  near  field,  but  the  compact  range  is  most  useful 
at  frequencies  above  1  GHz.  At  low  frequencies,  transverse  electromagnetic 
(TEM)  cells  [2]  and  parallel  plate  transmission  lines  are  used  to  produce  TEM 
fields  which  are  uniform  over  a  portion  of  their  cross  section.  Such  trans- 
mission line  structures  are  most  effective  for  testing  EUTs  that  are  elec- 
trically small.  There  exists  a  frequency  gap  (about  50  MHz  to  1  GHz)  between 
the  low  frequency  and  high  frequency  methods  where  another  method  of  producing 
a  plane  wave  is   needed. 

The  feasibility  of  using  phased  arrays  to  fill  the  frequency  gap  has  been 
studied,  and  a  near-field  array  synthesis  method  has  been  developed  to  compute 
the    array    element    excitations    [3].       The    method    optimizes    the    field    quality 
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throughout  the  test  volume  and  also  constrains  the  element  excitations  to 
reduce  the  field  strength  outside  the  test  volume.  The  synthesis  method  can 
be  applied  to  arbitrary  array  elements  and  geometries,  and  in  this  report  we 
apply  the  method  to  a  seven-element,  planar  array  of  dipole  antennas.  To 
check  the  theory  and  to  gain  experience  in  working  with  near-field  arrays,  we 
constructed  a  seven-element  array  of  dipoles  and  tested  it  at  a  frequency  of 
500  MHz.  In  an  actual  electromagnetic  susceptibility  test,  we  might  prefer 
array  elements  with  greater  bandwidth,  but  dipoles  were  useful  for  verifying 
the  array  theory  because  their  characteristics  are  well  understood.  Earlier 
experimental  results  have  been  obtained  by  Lynggaard  [4]  who  built  and  tested 
a  five-element  array  of  horns. 

The  organization  of  this  report  is  as  follows.  Section  2  contains  the 
analysis  and  numerical  results  for  the  seven-element  array.  Mutual  coupling 
between  the  array  elements  is  considered  in  section  2.4.  Section  3  describes 
the  array  construction  and  feed  network  and  the  method  of  probing  the  ampli- 
tude and  phase  of  the  electric  field.  Section  4  contains  the  experimental 
results  for  the  electric  field  measurements  throughout  the  test  volume  and 
gives  some  comparisons  with  the  theoretical  predictions  of  the  electric 
field.  Section  5  summarizes  the  results  of  this  study  and  makes  recommenda- 
tions  for  further  work. 

2.     Theory 

The  geometry  of  the  seven-element  planar  array  is  shown  in  figure  1.  The 
equally-spaced  array  elements  are  located  at  the  vertices  and  at  the  center  of 
a  regular  hexagon  in  the  x-z  plane.  This  configuration  can  also  be  viewed  as 
a  triangular  arrangement  of  the  array  elements,  and  it  is  known  to  be  very 
efficient  in  filling  the  planar  array  surface  without  producing  grating  lobes 
[5]. 

Only  the  driven  elements  are  shown  in  figure  1,  but  each  driven  element 
is  accompanied  by  a  passive  reflector  in  order  to  reduce  the  radiation  in  the 
backward  (negative  y)  direction.  Thus  each  array  element  is  actually  a  two- 
element,  Yagi-Uda  antenna  as  described  and  analyzed  in  section  2.1.  The  exci- 
tations of  the  Yagi-Uda  antennas  are  chosen  to  optimize  the  plane  wave  char- 
acter   of    the    field    in    the    rectangular    test    volume    shown    in    figure  1.       In 


section  2.2,  the  previously  developed  near-field  array  synthesis  method  [3]  is 
applied  to  this  array,  and  the  optimum  element  excitations  are  determined.  In 
section  2.3,  the  field  distributions   of  the  array  are  shown. 

2.1     Two-Element,   Yagi-Uda  Antenna 

Since  the  seven  array  elements  are  identical  and  mutual  coupling  is  ne- 
glected in  the  initial  analysis,  it  is  sufficient  to  analyze  a  single  isolated 
Yagi-Uda  antenna.  The  center  element  of  the  array  is  shown  in  figure  2,  and 
it  consists  of  a  driven  dipole  of  half  length  h^  and  a  reflector  of  half 
length  hp.  Both  dipoles  are  z  directed.  The  current  distributions  on  both 
the   driven   element    I^U')    ancl   the    reflector   element    Ir(z')    are  assumed  to   be 

sinusoidal : 

I  . 

I.(z')  =  -r-nTTT-  sin  kfh  .  -  Iz1  I] 

dv      '        sin  kh  .  ^   d       '      ' J 

a 

and  (1) 

I 

I   (z')   =      .    r, ,      sin   kfh     -   Iz'  II, 
rv      '       sin  kh  v    r       '      ' ; ' 

r 


where  k  =  w  /uoeo  =  2-n/X,  y0  and  e0  are  the  free  space  permeability  and  per- 
mittivity, and  X  is  the  free  space  wavelength.  The  time  dependence  exp  (jut) 
is  suppressed.  The  current  distributions  in  eq  (1)  are  symmetrical  in  z1  and 
are  zero  at  the  ends  of  the  elements.  A  more  precise  form  of  the  currents 
could  be  obtained  from  three-term  theory  [6,7]  or  the  method  of  moments  [8], 
but  the  sinusoidal  assumption  is  adequate  for  dipoles  of  length  on  the  order 
of  a  half  wavelength.  The  advantage  of  the  sinusoidal  current  assumption  is 
that  simple  expressions  are  available  for  the  near  fields  [9]  and  the  self  and 
mutual    impedances   [10]. 

The  base  currents   of  the  driven  and   reflector  elements,    1^  and   Ip,   can   be 
obtained  from  the  following  circuit   relations: 


V  ,  =  Z.I  ,  +  Z  ,   I 

d         d  d         dr  r 

and  (2) 

0  =  Z  ,   I  .  +  Z   I    . 
dr  d         r  r 


In  eq  (2),  Vd  is  the  voltage  applied  to  the  base  of  the  driven  element,  Zd  is 
the  self  impedance  of  the  driven  element,  Zp  is  the  self  impedance  of  the 
reflector  element,  and  Zdr  is  the  mutual  impedance  between  the  driven  and  re- 
flector elements.  The  voltage  applied  to  the  base  of  the  reflector  element  is 
zero  because  it  is  passive.  From  the  second  equation  in  eq  (2),  the  ratio  of 
the  reflector  current  to  the  driven  current  is 

Ir/Id  =  -Zdr/Zr.  (3) 

The  other  quantity  of  interest  is  the  input  impedance  of  the  driven  element 

zin 

V        Z  2 

S'n  "  I  .  Ld         Z  *  l4; 

d        r 

To  calculate  the  self  and  mutual  impedances,  we  use  the  mutual  impedance  ex- 
pressions given  by  King  [10].  These  expressions  involve  sine  and  cosine  inte- 
grals, and  efficient  computer  codes  are  available  for  evaluating  these 
integrals  [11].  Self  impedances  are  evaluated  by  computing  the  mutual  imped- 
ance at  a  separation  equal  to  the  dipole  radius  a. 

The  quantities  in  eqs  (2)  through  (4)  have  been  evaluated  for  the  fol- 
lowing parameters:  hd/X  =  0.25,  hr/x  =  0.26,  a/A  =  5.29  x  10"3,  and  d/X 
=  0.25.  The  value  of  a/X  was  determined  from  the  experimental  parameter  of  a 
=  0.318  cm  (0.125")  and  a  frequency  of  500  MHz.  For  the  above  parameters  the 
various  impedances  are: 

Zd  =  7.31  +  j  40.6  n, 

Zr  =  82.2  +  j  67.0  n,  and 

Zdr  =  43.2  -  j  29.8  a.  (5) 

By  substituting  eq  (5)  into  eqs  (3)  and  (4),  the  following  values  are  obtained 
for  the  current   ratio  and  the  input   impedance: 

Ir/Id  =  -  0.138  +  j  0.476 


and 

Zin  =  81.3  +  j   65.2  n.  (6) 

Thus  the  input  impedance  is  larger  than  the  self  impedance  of  the  driven  ele- 
ment given  in  eq  (5).  The  fields  of  the  two  elements  add  constructively  in 
the  forward  direction  (positive  y)  and  cancel  in  the  backward  direction  (nega- 
tive y)  because  the  phase  of  Ip/Ij  is  near  90  deg.  The  resultant  far-field 
radiation  patterns  in  the  two  principal  planes  are  shown  in  figure  3.  The  E 
plane  is  the  yz  plane,  and  the  H  plane  is  the  xy  plane.  The  patterns  were 
computed  from  the  far-field  expressions  derived  in  appendix  A.  The  front-to- 
back  field  ratio  is  2.94  (i.e.,  9.37  dB).  Some  calculations  were  done  for 
other  values   of  d  and   h_,   but  not  much   improvement  was   observed. 

In  the  following  sections,  the  near  fields  of  the  Yagi-Uda  antenna  are 
required.  For  the  sinusoidal  current  distribution  given  in  eq  (1),  the  elec- 
tric field   radiated  by  the  driven  element   E     is   given   by   [9] 


where 


rd  rd  rd  rd  ,.., 

E=xE+yE+zE,  (7) 

—  x       J     y  z    '  v    ; 


Ed  ,  x  d  =  y_  E 

x      p     p*    y      p     p' 


-jnIH  [z  -  h.J  expt-jkPx) 


[- 


p       4tt  sin  kh,   L       p  rx 

(z  +  h  ,)  exp(-jkr2)       2z  cos   kh  .  exp(-jkr0] 

+ _  9. 1 

P  r2  p  r0         J' 

d           -jnld          expt-jkrj       exp(-jkr2) 
Ez  =  4,  sin  khd   I Fi +  F2 "   2  cos  khd  exp(-jkr0)/r0], 


r,   =  /p2  +   (z  -  h    )2,     r2  =  /p2  +   (z  +  h    )2, 


rQ   =  /p2  +  z2,   p  =  /x2  +  y2,   and  r\  =  /y0/eQ. 
The  magnetic  field  radiated  by  the  driven  element   Hd  is   given  by   [9] 


Hd  =  x  Hd  +  y   Hd,  (8) 

x       J     y  v    ' 

where  Hd   =  -=£  1-T,       Hd  =  -  HT,  and 
x         p     y         y       p     <t>' 

jld 
H6  =  4ttP  sin  kh     [exPHkri)  +  exp(-jkr2)   -  2  cos   khd  exp(-jkr0)J. 

The  same  expressions  apply  for  the  fields  of  the  reflector  element,  _E_r  and  _H_r, 
if  we  replace  Id  by  Ir,  hd  by  hp,  and  y  by  y  +  d.  The  total  fields  of  the 
Yagi-Uda  antenna,  Ef  ahd  _H  ,  are  simply  the  sums  of  the  fields  of  the  driven 
and  reflector  elements. 

E t   =  Ed  +  Er  and  H*   =  Hd  +  Hr.  (9) 

The  fields  in  eqs  (7),  (8),  and  (9)  are  those  of  the  center  element  of  the 
array  in  figure  1,  but  the  expressions  are  also  applicable  to  any  of  the  six 
outer  elements  because  all  elements  are  identical.  It  is  only  necessary  to 
replace  x  by  x-X;  and  z  by  z-z^  where  the  element  is  centered  at  (x^z.-)  in 
the  xz  plane. 

2.2  Array  Element  Excitations 

In  this  section  we  use  the  previously  developed  near-field  synthesis 
method  [3]  to  determine  the  excitations  of  the  array  elements.  The  desired 
field  is  a  plane  wave  propagating  in  the  positive  y  direction: 

E=zE   andH=xH,  (10) 

-=p     pz    -f>     px' 

where  E   =  E0  exp  (-jky)  and  H   =  —  exp(-jky). 

The  array  and  test  volume  in  figure  1  have  two  planes  of  symmetry,  the  xy 
plane  and  the  yz  plane.  Consequently,  the  excitations  of  the  array  elements 
possess  the  same  symmetry  as  shown  in  figure  1,  and  there  are  only  three  in- 
dependent excitations.  I  is  the  base  current  of  the  driven  element  at  the 
origin,  I  is  the  base  current  of  the  four  driven  elements  above  and  below  the 
x  axis,  and  I  is  the  base  current  of  the  two  elements  on  the  x  axis  at  x  =  + 
s.  For  later  use  in  the  synthesis  procedure,  we  define  a  three-element  column 
matri  x: 


[I]  = 


II 
h 

I3 


(11: 


The  values  of  I  ,  I  ,  and  I  are  to  be  chosen  to  approximate  the  plane  wave 
given  by  eq  (10)  within  the  test  volume  shown  in  figure  1.  If  we  neglect 
mutual  coupling  between  the  array  elements,  then  the  corresponding  input  volt- 
ages (v  ,  v  ,  and  v  )  are  proportional  to  the  base  currents: 


(12: 


Vl 

Ii 

v2 

=  z.      • 
in 

I2 

_V3_ 

_I3_ 

where  Z^ n  is  given  by  eq  (4).   If  mutual  coupling  is  to  be  included,  then  Z- 
must  be  replaced  by  the  appropriate  impedance  matrix  as  shown  in  section 
2.4.   Of  course,  we  always  include  the  mutual  coupling  between  the  driven 
element  and  its  reflector  as  given  by  eqs  (2)  through  (4). 


The  first  step  in  the  synthesis  procedure  [3]  is  to  define  a  column 
matrix  [P]  of  the  desired  plane-wave  values  on  the  surface  S  of  the  test 
vol ume: 

Pi 


[P]  = 


P2 


^i 


(13) 


where 


P.  =  (E  .   -nn  x  H  )| 

— i   v-=^>  tan       -^p '  \r_  =  r. 


The  subscript  tan  indicates  the  component  tangential  to  S,  and  n  is  a  unit 
inward  normal.  The  element  _P.  is  a  tangential  vector  evaluated  at  the  point  _r_ 
=  _r\j  as  shown  in  figure  4.  M  is  the  number  of  points  considered  on  S.  In 
addition  we  define  an  M  x  3  rectangular  matrix  [T]: 


[T]  = 


Xl  1   Ii  2 


.13 


121   122   123 


Im1  Im2  Im3 


(14) 


where  Tj j   =   (E^  -  n  n  xH)^. 

The  fields  in  the  expression  for  JTjj  are  evaluated  for  I.-  =  1  and  all  other 
currents  equal  to  zero.  The  actual  expressions  for  the  electric  and  magnetic 
fields  of  the  Yagi-Uda  elements  are  given  in  eqs  (7)  through  (9).  The  reason 
for  using  the  hybrid  vector  involving  the  tangential  components  of  both  the 
electric  and  magnetic  fields  on  S  is  that  this  particular  hybrid  vector  pos- 
sesses the  desired  uniqueness  properties  when  specified  over  a  closed  surface 
[3]. 

The  synthesis  problem  amounts  to  determining  [I]  when  [P]  (the  desired 
plane-wave  field)  and  [T]  (the  solution  to  the  forward  problem)  are  known.  We 
can  write  the  approximate  relationship  between  [P]  and  [T][I]: 


[T][I]   -  [P]. 


(15) 


Since  [P]  is  an  M  x  1  column  matrix  of  two-component  (tangential)  vectors  and 
[I]  is  a  3  x  1  column  matrix  of  scalars,  eq  (15)  actually  represents  2M  equa- 
tions in  3  unknowns.  Normally,  2M  is  greater  than  3,  and  eq  (15)  can  only  be 
solved  approximately.  The  least  squares  solution  to  eq  (15)  yields  the  best 
approximation  to  the  plane-wave  field  on  the  surface  of  test  volume,  and 
uniqueness  guarantees  that  the  field  is  also  optimized  throughout  the  test 
volume   [3]. 

If  the  least  squares  solution  to  eq  (15)  yields  large  element  excitations 
and  large  fields  outside  the  test  volume,  a  more  useful  solution  is  obtained 
when  a  constraint  is  placed  on  the  source  norm  [3,12].  For  the  array  in 
figure  1,  it  is  natural  to  define  the  following  diagonal  weight  matrix  [A]  for 
the  source  inner  product: 


[A]  = 


*1 

0 

0 

0 

a2 

0 

0 

0 

as 

(16) 


where  a     =   1,   a     =  4,  and  a     =  2, 


The  values  of  the  elements  a,-  are  chosen  to  equal  the  number  of  array  elements 
with  base  current  I,-.  The  square  of  the  source  norm  |[I]||2  is  given  by 
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||[I]||2  =  [T]  [A][I]  =  I     a.|I.|2.  (17) 

j=l  J  J 

In  eq  (17),  *  denotes  complex  conjugate,  and  the  tilde  denotes  transpose.  The 
constraint  on  the  source  norm  is  written  [12] 

||[I]||2  <  C,  (18) 

where  C  is  a  positive  constant  to  be  chosen. 

The    least    squares    solution    to    eq   (15)    with    the    constraint    eq   (18)    re- 
quires the  eigenvalues  Xn  and  the  eigenvectors   [<j>]n  which  satisfy 

[T]*[T][*]n   =  Xn[A][4>]n.  (19) 

Matrix  multiplication  in  eq   (19)   has  the  usual   meaning,   and  the  dot   product   is 
understood  whenever  the  elements  multiplied  together  are  vectors.     The  matrix 
[T]   [T]  is   Hermitian,   and  efficient  computer  codes   [11]  are  available  for  com- 
puting   the    eigenvalues    and    eigenvectors.      The    orthogonality    relationship    of 
the  eigenvectors   i  s: 

Ri'wuij .  t  J; ; ;  •  •  <20) 

The  constrained  solution  for  [I]  is  [12]: 


n=l  n 


where  Cn  =  [?]*[T]*[P]. 


The  Lagrange  multiplier  a  satisfies 


c  -  1 


l,    JX     +  a) 
n=l  v  n    ' 


=  0. 


(22) 


The  left  side  of  eq  (22)  is  a  monotical  ly  increasing  function  in  a  which  is 
real  and  positive,  and  the  equation  is  easily  solved  numerically  for  a.  If  a 
is  set  equal  to  zero,  then  eq  (21)  gives  the  eigenvector  expansion  of  the 
usual  least  squares  solution  of  eq  (15). 

A  normalized  squared  error  En  can  be  defined  as 


where  [F]  =  [T][I]  = 


1  lf_-  -  M2 

|[F]  -  [Pill2  _  i=l  "         ' 
||[P]||Z    "  "M 

I    IP,- I2 
1=1  n 


Ll 


li 


^M 


(23) 


It  has  been  shown  [12]  that  as  C  is  decreased  both  a  and  En  increase.  How- 
ever, it  is  often  the  case  that  large  decreases  in  the  source  norm  j|I|i  can 
be  achieved  with  only  small  increases  in  En.  A  useful  value  of  C  can  be  ob- 
tained by  considering  a  uniform  array  where  1=1=1=  Iu.  In  this  case, 
C  is  given  by 


c-7|g*. 


(24) 


The  value  of  Iu  is  derived  in  appendix  B 


2.3  Numerical  Results 


A  computer  program  was  written  to  solve  for  [I]  by  eq  (21).  Because  of 
the  double  symmetry,  the  match  points  were  taken  over  only  one  quarter  of  the 
surface  S  of  the  test  volume  as  shown  in  figure  4.  The  number  of  points  in 
the  x,  y,  and  z  directions  are  Mx,  NL,  and  Mz,  respectively.   Consequently, 


the  spacings  of  the  points  in  the  x,  y,  and  z  directions  are  B/(2MY),  A/M 


r 
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and  C/(2MZ),  respectively.  The  numerical  results  were  found  to  be  insensitive 
to  Mx,  My,  and  Mz  for  spacings  less  than  or  equal  to  A/4.  There  are  MXMZ 
points  on  the  front  and  back  faces  (y  =  y  ±  C/2),  MyMz  points  on  the  side 
face  (x  =  B/2),  and  MXM  points  on  the  top  face  (z  =  A/2).  Thus  the  total 
number  of  points  M  is 

M  =  2M  M  +  M  (M  +  M  ).  (25) 

x  z    yx  x    z'  v  ' 

The  numerical  results  in  figures  5  through  15  apply  to  the  following 
array  parameters:  s  =  A,  hd  =  d  =  A/4,  hp  =  0.26  A,  and  a  =  5.29  x  10"3  A. 
In  figures  5  through  11,  the  test  volume  is  a  1.5  A  cube  (A  =  B  =  C  =  1.5  A) 
which  is  centered  1.75  A  from  the  array  (y  =  1.75  A).  Thus  the  front  of  the 
test  volume  is  one  wavelength  from  the  array  (y  -  C/2  =  A).  The  magnitudes 
and  phases  of  the  synthesized  currents  are  shown  in  figure  5.  For  this  par- 
ticular case,  the  constraint  has  no  effect  because  the  unconstrained  currents 
(a  =  0)  satisfy  eq  (18).  The  currents  in  figure  5  are  normalized  to  the  cur- 
rent of  the  uniform  array  I  as  given  in  appendix  B.  The  center  current  I 
has  the  largest  magnitude,   and  this  was  true  for  all    cases  which  were  studied. 

For  the  synthesized  currents  in  figure  5,  the  normalized  squared  error  E 
as  given  by  eq  (23)  is  0.077  compared  to  En  =  0.385  for  the  uniform  array. 
This  difference  in  the  quality  of  both  the  magnitude  and  phase  of  the  electric 
field  is  evident  in  figures  6  through  11.  In  all  cases  the  plane-wave  value 
is  shown  over  the  extent  of  the  test  volume  for  comparison.  The  decay  in  the 
field  for  x  and  z  outside  the  test  zone  in  figures  7  and  8  is  about  the  same 
for  the  synthesized  and  uniform  arrays.  In  the  z  =  0  plane,  the  electric 
field  has  only  a  z-component,  and  there  are  no  cross-polarized  components. 
For  z  *  0,  there  is  a  small  cross-polarized  field  (Ey)  as  shown  in 
figure  12.  (Here  we  refer  to  either  orthogonal  component  of  the  electric 
field,  Ex  or  E  ,  as  a  cross-polarized  component.  In  propagation  theory,  E  is 
more  often  called  a  longitudinal  component  producing  a  wave  tilt.)  In  general 
the  amplitude  errors  in  figures  6  through  8  tend  to  be  more  significant  than 
the  phase  errors  in  figures  9  through  11.  Although  not  shown  here,  the  errors 
in  the  magnetic  field  have  similar  characteristics. 

There   are   several    ways   to   reduce  the   size  of  the  field  errors:      (1)   in- 
crease the  number  of  array  elements,   (2)  decrease  the  size  of  the  test  volume, 
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(3)  move  the  test  volume  farther  from  the  array.  An  increase  in  the  number  of 
elements  gives  more  degrees  of  freedom  for  the  synthesis  procedure,  but 
increases  the  cost  and  complexity  of  the  array.  A  decrease  in  the  size  of  the 
test  volume  makes  it  easier  to  maintain  a  uniform  field,  but  limits  the  size 
of  the  test  object.  In  figures  13  through  17,  we  show  numerical  results  for 
y  =  2.75  X.  Thus  the  front  face  of  the  test  volume  is  now  two  wavelengths 
from  the  array.  All  other  parameters  remain  the  same.  The  synthesized  cur- 
rents both  with  and  without  a  constraint  are  shown  in  figure  13.  For  this 
case,  the  normalized  squared  error  En  is  0.020  for  the  unconstrained  synthe- 
sized array,  0.046  for  the  constrained  synthesized  array,  and  0.340  for  the 
uniform  array.  Thus  the  errors  are  smaller,  as  expected,  when  the  test  volume 
is  farther  from  the  array.  The  magnitude  of  the  electric  field  is  shown  along 
the  three  axes  in  figures  14  through  16,  and  the  errors  are  seen  to  be  smaller 
than  those  in  figures  6  through  8.  The  cross-polarized  component  is  shown  in 
figure  17,  and  it  is  seen  to  be  smaller  than  in  figure  12. 

Because  it  is  difficult  to  display  fields  throughout  a  volume,  only  the 
fields  on  the  three  centerlines  of  the  cube  have  been  shown.  However,  the 
normalized  squared  error  En  which  represents  the  surface  field  errors  is  a 
good  indicator  of  average  field  errors  throughout  the  volume. 

Some  cases  were  run  for  different  values  of  the  element  spacing  as  shown 
in  figure  1,  and  the  best  results  were  obtained  for  s  =  X.  For  smaller  values 
of  s,  the  array  size  gets  smaller,  and  it  becomes  necessary  to  decrease  the 
size  of  the  test  volume.  For  larger  values  of  s,  undesirable  grating  lobes 
begin  to  form.  If  we  assume  that  spacings  on  the  order  of  one  wavelength  are 
good  for  broadside  arrays  which  are  not  electronically  scanned,  then  the  width 
of  the  test  volume  must  be  somewhat  less  than  /N  X,  where  N  is  the  number  of 
elements.  In  other  words,  the  width  of  the  array  must  be  somewhat  greater 
than  the  width  of  the  test  volume. 

2.4  Mutual  Coupling 

In  the  previous  sections,  the  optimum  values  of  the  base  currents  were 
determined.  When  mutual  coupling  between  the  array  elements  is  neglected,  the 
input  voltages  are  proportional  to  the  base  currents  as  shown  in  eq  (12).  In 
this  section  we  derive  more  accurate  expressions  for  the  input  voltages  which 
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include  the  effect  of  mutual  coupling.   The  results  tell  us  when  mutual 
coupling  is  important  and  how  to  account  for  it. 

The  base  currents  I,-  and  input  voltages  v^  for  the  seven-element  array 
are  shown  in  figure  18.  The  voltages  and  currents  satisfy  the  following 
matrix  relationship: 


Vl 

v2 

• 

= 

• 

• 

v7 

Zh 

Z21 
Z71 


Z12 

Z22 


17 


•77 


(26) 


The  mutual  impedances  Z- •  are  defined  as 

v. 


(27) 


where  all  base  currents  except  I-  are  zero 


The  reflector  currents  are 
assumed  to  satisfy  the  approximate  relationship  given  in  eq  (3).  Thus  the 
mutual  impedance  in  eq  (27)  can  be  written  as  the  sum  of  two  terms: 


I 


z..  =  1.: 


^zi 


(28) 


where  Z^  •  is  the  mutual  impedance  between  the  i  and  j  driven  elements  and 
Z^^  is  the  mutual  impedance  between  the  i  driven  element  and  the  j  re- 
flector. The  ratio  (I  /I^)  is  approximately  given  by  eq  (3).  A  more  accurate 
expression  for  (Ir/I^)  can  be  obtained  by  adding  seven  more  equations  to  eq 
(26)  and  solving  the  entire  set  of  equations  simultaneously  [7],  but  we  have 
found  that  the  resultant  value  for  (I./Ij)  varies  little  from  that  given  by 
eq  (3).  The  expressions  for  Z^  _•  and  Zn-  -r  in  eq  (28)  are  those  for  dipoles  in 
echelon  [10], 

Because  of  the  double  symmetry  of  the  array,  many  of  the  voltages  and 
currents  in  figure  18  are  equal: 

v2  =  vH  =  v5  =  v7,  v3  =  v6, 
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I2  =  I4  =  I5  =  I7»  and  I3  =  I6. 


(29) 


Thus  there  are  only  three  independent  voltages  (v  ,  v  ,  v  )  and  currents  (I 
I  ,  I  )  as  shown  ii 
duced  to  three  rows 


I  ,  I  )  as  shown  in  figure  1,  and  the  matrix  equation  in  eq  (26)  can  be  re- 
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Many  of  mutual  impedances  are  equal : 


and 


zn   z22  z33, 

7   =7   =7   =7   =7   =7   =7   =Z 
^12    ^14    Llb         L\7         ^21    ^"23    ^"32    ^i*' 


Z13  ::  Z16    Z27  ~  Z3i  . 


(30) 


Z26  =:  Z37  "  Z35* 


(31) 


By  substituting  eq  (31)  in  eq  (30),  we  have 


Vl 

v2 
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Zn  4  Z12        2  Z13 
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Z  in  eq  (38)  is  equal  to  Zq- n  in  eqs  (4)  and  (12).  Thus,  eq  (32)  reduces  to 
eq  (12)  only  when  the  mutual  impedances  U  »  Zi3,  Z  ,  Z25>  Z26,  and  Z3fe)  are 
all  zero.  These  mutual  impedances  decrease  in  magnitude  as  the  interelement 
spacing  s  is  increased. 

Since  the  currents  (I  ,  I  ,  I  )  are  already  known  from  eq  (21),  the  re- 
quired input  voltages  (v  ,  v  ,  v  )  can  be  obtained  directly  from  eq  (32)  by 
matrix  multiplication.  The  resultant  input  voltages  are  shown  in  figures  19 
and  20  for  the  same  parameters  which  were  used  to  obtain  the  currents  in 
figure  5  and  the  constrained  currents  in  figure  13.  The  voltages  are  normal- 
ized to  the  voltage  vu  for  a  uniform  array  which  is  given  by 

vu  =  Iu  Zin,  (33) 

where  Iu  is  given  by  eq  (B-7)  and  Z.,-n  is  given  by  eq  (4).  The  voltages  in 
figures  19  and  20  which  do  not  include  mutual  coupling  were  computed  from  eq 
(12).  Thus  the  differences  between  the  dashed  and  solid  results  are  an  indi- 
cation of  the  importance  of  mutual  coupling.  The  differences  occur  primarily 
in  the  phases,  but  this  is  not  necessarily  true  for  other  cases.  The  effect 
of  the  interelement  spacing  s  on  mutual  coupling  has  been  studied,  and  the 
mutual  coupling  has  been  found  to  decrease  as  s  is  increased  as  expected. 
This  is  why  the  hexagonal  geometry  was  chosen  to  allow  a  larger  interelement 
spacing  [5], 

3.  Experiment 

The  array  as  described  in  section  2  and  shown  in  figures  1  and  2  was  con- 
structed, and  a  series  of  measurements  was  performed  throughout  the  test 
volume  shown  in  figure  1.  A  test  frequency  of  500  MHz  was  chosen  for  conven- 
ience, although  the  array  can  be  adjusted  for  any  frequency  from  100  to 
1000  MHz. 

3.1  Array  Construction 

Each  of  the  seven  two-element  Yagi-Uda  antennas  was  constructed  as  shown 
in  figure  21.  The  driven  dipole  element  is  fed  with  a  dual  coaxial  trans- 
mission line  as  described  in  [13],  The  coaxial  lines  were  in  turn  fed  with  a 
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0°  -  180°  hybrid  junction  (power  splitter),  as  shown  in  figure  22,  to  provide 
optimum  excitation  of  the  dipole.  Each  dipole  antenna  was  then  tuned  for 
resonance  at  500  MHz  by  adjusting  the  length  of  the  telescoping  arms  on  the 
dipole.  The  driven  to  reflector  distance  d  and  the  reflector  length  2hp 
(figure  2)  were  then  adjusted  according  to  dimensions  in  section  2.3;  d  =  A/4 
and  hp  =  0.26  X,  based  on  the  physical  length  of  the  driven  element  (X/2). 

The  masts  for  the  seven  Yagi-Uda  antennas  were  attached  to  the  structure 
shown  in  figure  23.  This  particular  arrangement  allowed  easy  positioning  of 
all  the  elements  for  any  chosen  frequency  (100  to  1000  MHz),  with  500  MHz 
chosen  for  this  experiment.  The  dual -coaxial  lines  from  each  antenna  extended 
back  and  perpendicular  from  the  driven  element  to  a  wall  of  rf  absorber 
(figure  24).  The  rf  source,  power  splitting  network,  and  instrumentation  are 
behind  the  absorber  wall  and  outside  the  anechoic  chamber. 

The  feed  system  (power  splitting  network)  consisted  of  13  broadband 
hybrid  junctions,  two  phase  shifters,  and  two  variable  attenuators.  Six  junc- 
tions were  connected  as  in-phase  power  splitters  and  seven  were  used  to  feed 
the  dual -coaxial  lines  from  the  antennas;  hence  these  were  connected  to 
utilize  the  0°  and  180°  ports.  The  system  is  better  described  in  figure  25. 

3.2  Probe 

The  rf  probe  construction  duplicated  the  driven  element  of  the  Yagi-Uda 
antennas.  Here  the  dipole  length  was  1/10  X  (approximately  6  cm)  rather  than 
X/2  used  for  the  resonant  antenna.  The  same  dual-coaxial  feed  system  was  used 
and  these  cables  were  carefully  positioned  to  minimize  interactions  with  the 
rf  field  being  measured.  As  near  as  was  practical,  the  cables  were  hidden 
behind  sections  of  rf  absorber.  The  signal  from  the  probe  then  traveled 
through  coax  line  along  the  floor  of  the  anechoic  chamber  and  out  through  the 
door  to  the  instrumentation.  Both  amplitude  and  phase  of  the  electric  field 
components  were  measured. 

3.3  Track 

The  probe  moved  through  the  test  volume  on  a  nonconducti ve  motorized 
track.  The  motor  for  the  track  sat  on  the  base  and,  through  a  series  of  belts 
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and  pulleys,  moved  the  probe  in  a  horizontal  linear  direction.  The  motor  also 
drove  a  precision  potentiometer.  The  resistance  of  the  potentiometer  was  then 
used,  after  careful  calibration,  to  determine  the  position  of  the  probe  along 
the  track.  The  details  of  the  track  may  be  seen  in  the  photograph  in  figure 
26.  The  direction  of  the  track  was  aligned  with  the  z  axis  of  the  array 
coordinates,  thus  all  scans  are  along  this  axis.  The  x  position  can  be 
changed  by  moving  the  cart  (a  feature  of  the  anechoic  chamber)  along  the  y 
coordinate,  perpendicular  to  the  array.  Thus  by  selecting  an  x  and  y  position 
and  scanning  the  z  direction,  the  test  volume  can  be  mapped.  The  amplitude 
and  phase  measurements  were  taken  by  manually  controlling  the  track  position 
via  a  switch  located  outside  the  anechoic  chamber.  At  every  selected  loca- 
tion, the  computer  was  triggered  to  accumulate  the  data  and  store  it  on 
disk.  A  single  scan  read  60  to  80  data  points.  A  more  elaborate  arrangement 
would  have  allowed  the  computer  to  control  the  track  also;  however,  this  sys- 
tem proved  adequate  for  these  measurements. 

3.4  Instrumentation 

The  measurement  system  (fig.  24),  used  a  vector  analyzer  (phase  amplitude 
receiver),  digital  multimeter  to  read  ohms  and  hence  track  position,  and  a 
computer/instrument  controller.  As  indicated,  the  track  was  positioned  manu- 
ally but  all  the  information  was  read  with  the  computer.  The  computer  read 
both  the  field  parameters  and  the  track  position  at  every  test  point.  The 
data  were  then  saved  on  floppy  disk  for  later  processing. 

4.  Experimental  Results  and  Comparisons  with  Theory 

Four  sets  of  element  weightings  were  applied  to  the  seven  element 
array.  These  amplitude  and  phase  feeding  schemes  were  calculated  using  the 
synthesis  techniques  outlined  in  section  2.  The  calculations  used  two  cubic 
(A  =  B  =  C  =  1.75  X  in  figure  1)  test  volumes  at  different  distances  from  the 
array,  each  with  and  without  mutual  coupling  adjustments.  We  identify  the 
four  sets  of  element  weightings  as  feeds  1-4  described  below: 

Feed  1.  Test  volume  centered  at  x  =  0,  y  =  1.75  X,  and  z  =  0  without  correc- 
tions for  mutual  coupling  between  array  elements  (fig.  19). 
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Feed  2.  Same  as  1,  but  with  corrections  for  mutual  coupling  between  array 
elements  (fig.  19). 

Feed  3.  Test  volume  centered  at  x  =  0,  y  =  2.75  x,  and  z  =  0  without  correc- 
tions for  mutual  coupling  between  array  elements  (fig.  20). 

Feed  4.  Same  as  3,  but  with  corrections  for  mutual  coupling  between  array 
elements  (fig.  20). 

The  test  results  are  shown  in  figures  27  to  53.  In  all  cases,  the  mea- 
sured quantity  is  shown  as  a  function  of  the  z  coordinate,  f(z).  Figures  27 
through  41  show  both  theoretical  and  measured  values.  Figures  42  to  53  are 
measured  values  only  and  show  more  detail  throughout  the  test  volumes.  All 
figures  are  scans  parallel  to  the  z  axis  at  some  distance  y  from  the  array  and 
lying  in  an  x  =  constant  plane.  Feed  1  was  probed  in  the  x  =  -0.75  X,  x  =  0, 
and  x  =  +0.75  X  planes,  while  feeds  2,  3,  and  4  were  probed  only  in  the  x  =  0 
plane.  All  amplitude  and  phase  values  are  referenced  to  the  center  of  the 
test  volume;  values  there  are  0  dB  and  0  deg.  The  test  volume  boundaries  are 
shown  by  the  vertical  columns  of  dots. 

Figures  27  through  29,  respectively,  show  scans  in  the  plane  x  =  0  along 
the  near  face  (y  =  1.0  X),  center  (y  =  1.75  X),  and  far  edge  (y  =  2.5  X)  of 
the  test  volume  used  for  feed  1.  All  three  electric  field  components  (Ex,  E  , 
and  Ez)  were  calculated  and  measured  for  each  of  these  locations.  The  calcu- 
lated curves  (shown  by  the  symbol  +)  are  obviously  not  those  of  a  perfect 
plane  wave;  rather  they  represent  the  best  approximation  to  the  plane  wave  as 
determined  by  the  synthesis  technique  for  an  array  of  only  seven  elements.  As 
the  number  of  elements  increases,  the  field  becomes  closer  to  that  of  a  plane 
wave  [3]. 

A  uniform  plane  wave  would  have  a  constant  Ez  with  no  Ex  or  E  compo- 
nents; hence,  Ex  and  E  may  be  considered  error  fields  for  this  experiment. 
Theory  predicts  an  Ey  component  at  these  close  distances  to  the  array  and  mea- 
surement confirms  it.  E  is  strongest  closest  to  the  array  (y  =  1.0  X)  and 
decreases  as  the  distance  increases.  Ex  is  predicted  to  be  zero;  the  measure- 
ment shows  that  it  is  small,  but  not  zero  (10  to  20  dB  or  more  below  Ey). 
These  data  may  indicate  the  effectiveness  of  the  rf  probe  in  being  able  to 
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resolve  the  field  components  or  other  measurement  errors.  Hence,  little  or  no 
emphasis  should  be  placed  on  interpretation  of  the  Ex  data. 

The  Ex  and  E  components  are  low  enough  as  to  have  an  insignificant 
effect  on  typical  EM  susceptibility  measurements.  So  the  remaining  figures 
(30  to  53)  detail  only  the  z-di  rected  electric  field  Ez.  The  figures  are 
paired  to  show  amplitude  and  phase  under  the  same  conditions;  i.e.,  figures  30 
and  31  show  amplitude  and  phase  for  feed  1  in  the  x  =  0  plane  at  y  =  1.0  X, 
1.75  X,  and  2.5  X,  figures  32  and  33  show  the  same  for  the  x  =  -0.75  X  plane, 
and  figures  34  and  35  are  for  the  x  =  +0.75  X  plane.  The  amplitude  data,  as 
presented,  can  be  interpreted  as  follows.  Each  grid  line  directly  above  any 
given  pair  (calculated  and  measured)  of  curves  for  a  particular  y  value  cor- 
responds to  0  dB.  For  example,  figure  30  shows  three  pairs  of  curves,  one 
each  for  y  =  1.0  X,  y  =  1.75  X,  and  y  =  2.5  X.  The  center  value  for 
y  =  1.75  X  is  0  dB  as  this  is  the  reference  point.  The  center  value  for 
y  =  2.5  X  is  near  0  dB  for  the  measured  and  about  -0.7  dB  for  the  calculated 
value.  The  values  for  y  =  1.0  x  at  the  center  are  about  -3  to  -4  dB.  Again, 
these  values  are  measured  from  the  grid  line  directly  above  the  curves,  with 
10  dB  between  grid  lines.  It  is  interesting  to  note  that  the  field  amplitude 
does  not  fall  off  as  1/r  as  one  proceeds  away  from  the  array,  but  rather 
increases  initially  and  then  decreases  more  slowly  than  1/r.  This  trend  is 
shown  in  more  detail  for  the  theoretical  results  in  figures  6  and  14.  Feed  2 
results  are  shown  in  figures  36  and  37  for  the  x  =  0  plane. 

Feed  3  refers  to  the  more  distant  test  volume,  and  figures  38  and  39  show 
data  for  the  x  =  0  plane  at  the  near  face  (y  =  2.0  X),  center  (y  =  2.75  X), 
and  the  far  boundary  (y  =  3.5  X)  of  the  test  volume.  The  final  pair  of  graphs 
with  both  theoretical  and  measured  values  compared  (figs.  40  and  41)  detail 
feed  4  for  the  same  locations  as  feed  3. 

The  agreement  shown  between  experiment  and  theory  in  figures  27  to  41  is 
quite  good.  There  is  one  set  of  data  (figs.  32  and  33),  howver,  where  an 
apparent  error  in  positioning  the  probe  caused  a  shift  in  the  experimental 
data  for  y  =  1.0  X.  Feeds  2  and  4  were  used  to  account  for  mutual  coupling 
between  the  elements,  but  mutual  coupling  effects  are  so  small  that  little 
improvement  is  seen  in  figures  36,  37,  40,  and  41.  Mutual  coupling  would  be 
more  important  if  the  elements  were  more  closely  spaced. 
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In  general  the  measurement  agrees  well  enough  with  theory  that  the  re- 
maining measurement  data  for  additional  locations  are  presented  alone. 
Figures  42  and  53  repeat  the  same  locations  as  shown  in  figures  27  to  41  with 
additional  traces  added  to  show  every  1/4  wavelength  in  y  for  each  chosen  x 
plane.  Again,  the  amplitude  and  phase  values  are  referenced  to  the  center  of 
the  test  volume,  the  amplitude  traces  have  a  0  dB  line  directly  above  each 
one,  and  the  corresponding  phase  for  each  amplitude  trace  is  on  the  following 
figure.  This  density  of  traces  is  shown  to  indicate  the  field  uniformity 
along  the  y  dimension  in  addition  to  the  z  scan.  The  amplitude  curves  can  be 
visually  grouped  to  show  a  third  dimension  and  give  the  viewer  a  feel  for  the 
field  intensity  anywhere  in  the  constant  x  planes. 

5.  Conclusions  and  Recommendations 

A  general  near-field  array  synthesis  technique  [3]  has  been  applied  to  a 
seven-element  array  of  Yagi -Uda  antennas  in  order  to  optimize  the  field  uni- 
formity throughout  a  rectangular  test  volume.  The  array  was  constructed  and 
tested  in  an  anechoic  chamber  at  a  frequency  of  500  MHz  to  verify  the  theory 
and  to  gain  experience  in  controlling  the  amplitude  and  phase  of  the  array 
elements.  The  agreement  between  the  measurements  and  theory  was  excellent  for 
both  the  amplitude  and  phase  of  the  electric  field.  The  effect  of  mutual 
coupling  between  the  array  elements  was  sufficiently  small  that  the  feed  cor- 
rections to  account  for  mutual  coupling  produced  a  negligible  change  in  the 
measured  fields. 

There  are  a  number  of  logical  extensions  to  this  work.  The  magnetic 
field  could  be  probed  throughout  the  test  volume  with  a  small  loop  probe.  The 
Yagi-Uda  array  was  chosen  for  convenience  of  analysis  and  construction  in 
validating  the  theory.  However,  to  make  the  phased  array  useful  for  electro- 
magnetic susceptibility  testing,  we  need  to  be  able  to  scan  the  frequency  over 
a  large  bandwidth,  perhaps  one  or  two  octaves.  This  probably  can  be  achieved 
by  computing  frequency-dependent  weightings  for  the  array  elements,  but  di- 
poles  might  not  be  the  most  useful  elements  because  of  their  frequency- 
dependent  input  impedance.  To  achieve  very  wide  bandwidths  on  the  order  of  a 
decade,  wide-band  elements  would  probably  be  required.  In  order  to  avoid 
having  to  rotate  the  test  object,  the  angle  of  arrival  of  the  plane  wave  could 
be  scanned  electronically,  and  angular  scanning  would  require  a  larger  number 
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of  elements.   Also,  the  appropriate  geometry  for  angular  scanning  would  re- 
quire further  study. 
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Appendix  A.     Far-Field  Pattern  of  a  Two-Element,   Yagi-Uda  Antenna 

In  the  far  field,  the  Yagi-Uda  antenna  of  figure  1  radiates  only  a  theta 
component  of  electric  field.  In  spherical  coordinates  (r,  e,  <(>),  the  far- 
field  of  the  driven  element  EQd  is  given   by  [9] 

,       jnl  .  exp(-jkr)     cos(kh  ,  cos  e)   -  cos  kh  , 

F  d  = 9 [ 9 9.1  (i\_i) 

LQ  2ttt  sin  kh,       L  sine  J*  lM  l) 

a 


The  far-field  of  the  reflector  element  EQr  is  given  by 


jnl  .  exp(-jkr  )  cos(kh  cos  e)  -  cos  kh 

E   r  =  9 L_  f C Hi  ri_o) 

Le  2ttt  sin  kh  L  sin  e  J'  K       ' 

r 


where  rr  =  r  +  d  sin  8  sin  <f>.     The  total    far-field  Efit  can  be  written 


t       r  d       r  r       jnId  GXP("J'kr) 


Ee     "  V  +  Ee     "  2,r P<6'*>  ^ 


where 


,         cos(kh.cose)  -  cos  kh . 


sinkh  , 
d 


I     [cos(kh  cose)   -  cos  kh   ]exp(-jkd  sine  sin<j>) 

+  s. C C 1 

I.  sinkh  '* 

d  r 


The  total   magnetic  field  H.^  is  given  by 


H  *  =  EgVn.  (A-4) 

The  pattern  in  the  E  plane  (yz  plane)  is  obtained  from  P(e,  tt/2),  and  the 
pattern  in  the  H  plane  (xy  plane)  is  obtained  from  P(tt/2,  <f>).  For  the  case 
plotted  in  figure  3  where  d  =  x/4  and  the  phase  (lr/l(i)  is  near  tt/2,  the  maxi- 
mum of    |P|   occurs  at  e  =  <j>  =  it/2. 
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Appendix  B.  Uniform  Array 

In  order  to  determine  the  actual  value  of  the  base  current  for  uniform 
excitation  Iu,  we  first  consider  the  equivalent  smoothed  current  density^. 
From  the  previous  results  [3],  J,,  is 


2  E 


J  =  2  y  x  H  I  n 
-s    J      -p|y=0 


o  - 
-  z. 


(B-l) 


The  effective  dipole  moment  of  each  element  ^_  shoul d  then  equal  the  sheet  cur- 
rent times  the  area  of  a  hexagonal  cell,  /3  s2/2: 


v^  s2  E0  . 
£= « — z- 


(B-2) 


The  effective  dipole  moment  of  the  driven  element  jd^  is  given  by 


^  -z    j     I  <z')  dz', 

-hd 


where 


(B-3) 


I.(z')  =  I 
dv  '    u 


sin  k(hd  -  |z|) 
sin  kh. 


The  integration  in  eq  (B-3)  can  be  performed  to  yield 

2(1  -  cos  kh.)  „ 
^d  ''    lu         k  sin  kh,   z* 


(B-4) 


Similarly,  the  effective  dipole  moment  of  the  reflector  element  _p_p  is 

2(1  -  cos  khr)  A 
£-r  =  1r         k  sin  kh    Zo 


(B-5) 


where  Ip  =  -Iu  l^r/lr.  Since  the  reflector  parameters  are  chosen  to  yield  an 
in-phase  contribution  in  the  forward  direction,  the  total  effective  dipole 
moment  jjj.   i  s 


fit 


2   1       (1  -  cos  kh  ,) 

U      r  V 


sin  kh 


dr 


(1   -  cos  kh   ) 

LLl 

sin  kh         J* 


(B-6) 
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Iu  is  determined  by  equating  jdj.   in  eq   (B-6)  with  g_  in  eq   (B-2) 


I     = 
u 


-E0ks2/3     (1   -  cos   khj 


2n 


sin  kh 


"dr 


(1  -  cos  kh   )     n 

Lli-1 

sin  kh         J 


(B-7) 
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Figure  1.        Seven-element    (equally   spaced)    array   in  the   xz   plane.     The 

rectangular  test   volume   (A  x  B  x  C)   is  centered  on  the  y-axis  at 
=  y0._    Because  of  symmetry,  there  are  only  three  independent 


excitations    (1,1 


L2' 


V- 
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Figure  2.       Two-element,    Yagi-Uda  antenna.     Both  the  driven  and   reflector 
elements   are  cylindrical    dipoles   of   radius   a. 
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Figure  3. 


Normalized 
antenna  in  the  two 
h  /A  =  0.26,   d/X   = 


far-field  pattern 
principal 
0.25,   and 


of  a  two-element,   Yagi-Uda  antenna 


planes.     Parameters: 
a/X  =   5.29  x   10-3 


hd/X  =  0.25, 
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Figure  4. 


Closed  surface  S  of  test   volume  with. two  planes   of  symmetry.     The 
match  points  are  located  at   r, ,   and  n  is  the  unit  inward  normal. 
For  this   case,   M     =  M     =   2  aTTd  M     =  4. 
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Figure  5.   Magnitude  and  phase  of  the  synthesized  currents 
the  constraint  has  no  effect. 
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Figure  6.   Normalized  magnitude  of  the  electric  field  for  a  synthesized  array 
and  a  uniform  array.  The  perfect  plane  wave  is  shown  for 
comparison. 
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Figure  7.   Normalized  magnitude  of  the  electric  field  for  a  synthesized  array 
and  a  uniform  array.  The  perfect  plane  wave  is  shown  for 
comparison. 
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Figure  8.   Normalized  magnitude  of  the  electric  field  for  a  synthesized  array 
and  a  uniform  array.  The  perfect  plane  wave  is  shown  for 
comparison. 
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Figure  9.   Phase  of  the  electric  field  for  a  synthesized  array  and  a  uniform 
array.  The  perfect  plane  wave  is  shown  for  comparison. 
Parameters:  A  =  B  =  C  =  0,  yQ  =  1.75X,  and  x  =  z  =  0. 
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Figure  10.  Phase  of  the  electric  field  for  a  synthesized  array  and  a  uniform 
array.  The  perfect  plane  wave  is  shown  for  comparison. 
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Figure  11.  Phase  of  the  electric  field  for  a  synthesized  array  and  a  uniform 
array.  The  perfect  plane  wave  is  shown  for  comparison. 
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Figure  12.     Normalized  cross-polarized  component  of  electric  field  for  a 
synthesized  array  and  a  uniform  array. 
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Magnitude  and  phase  of  the  synthesized  currents  with  and  without  a 
constrai  nt. 
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Figure  14.  Normalized  magnitude  of  the  electric  field  for  a  synthesized 
array,  a  constrained  synthesized  array,  and  a  uniform  array, 
perfect  plane  wave  is  shown  for  comparison. 
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Figure  15.  Normalized  magnitude  of  the  electric  field  for  a  synthesized 

array,  a  constrained  synthesized  array,  and  a  uniform  array.  The 
perfect  plane  wave  is  shown  for  comparison. 
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Figure  16.  Normalized  magnitude  of  the  electric  field  for  a  synthesized 
array,  a  constrained  synthesized  array,  and  a  uniform  array, 
perfect  plane  wave  is  shown  for  comparison. 
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Figure  17.     Normalized  cross-polarized  electric  field  for  a  synthesized  array, 
a  constrained  synthesized  array,   and  a  uniform  array. 
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Figure  18.     Voltages   and  currents   for  the  seven-element   array. 
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Figure  20.  Magnitude  and  phase  of  the  voltages  with  and  without  mutual 

coupling.  These  voltages  correspond  to  the  currents  in  figure  13, 
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Figure  23.     Seven-element   array  and  support   structure, 


47 


RF 
SOURCE 


FEED 
SYSTEM 


REFERENCE 
SIGNAL 


TRRCK 
CONTROL 


LEFT 


RIGHT 


RNECHOIC  CHAMBER 


7  ELEMENT  ORRRY 
e — 0 


-e — O 


-e— o 


-e— O 


-e— e 


COMPUTER 


RF    PROBE 


TRACK    — 


CRRT 


O      O 


DIGITAL 
MULTIMETER 


VECTOR 
ANALYZER 


Figure  24.  Measurement  system. 


48 


-t-> 
in 
>> 
(/> 

■a 
ai 
a; 


<T3 


CM 


ai 


49 


Figure  26.     Experiment  hardware: 
seven-element  array. 


probe  positioner   (track  mechanism)   and 


50 


RMPLITUDE  (Feed  #1) 


10 


0 


-10 


m 


u 

Q 

D 
h 

M 

J 

Q_ 

E 
d 


-20 


-30 


-40 


-50 


60 


-2X 


1 r—    i 1 1 

i — i — i — 

1 

-T      n 1 ] 1 1 1 

.   Theo  ry     +  +  +  +  + 
,  Me  asu  red 

^♦-vEz^^K 

I  / 

II           \          II          \          r                  \\ 

II 

\  Vy       Ex^ 

\ 

f'         \  II  ■            H 

i          Wi  •              \ 

\ 

r\  r\  A-l 

\  /  '      ^     \l 

.    .    ,    i 

1 1 1 

1 

i    i    i    i    i    j..  _i 

— X  Center  X 

Scan     the     Z     dimension 

Data     -For     f  CZ)  ,      X  =  0  ,      Y=  1  .  0 


2X 


Figure  27.     Theoretical    and  measured  electric  field  amplitude  at  the  center 
surface  of  the  test   volume    (y  =  X). 
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Figure  28.     Theoretical    and  measured  electric  field  amplitude  at  the  far 
surface  of  the  test   volume   (y  =   1.75  X). 
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Figure  29.     Theoretical    and  measured  electric  field  amplitude  at  the  near 
surface  of  the  test   volume   (y  =   2.5  A). 
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Figure  30.  Theoretical  and  measured  amplitude  of  Ez  at  three  distances  from 
the  array. 
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Figure  31.  Theoretical  and  measured  phase  of  Ez  at  three  distances  from  the 
array. 
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Figure  32.  Theoretical  and  measured  amplitude  of  Ez  at  three  distances  from 
the  array  for  x  =  0.75  X. 
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Figure  33.  Theoretical  and  measured  phase  of  Ez  at  three  distances  from  the 
array  for  x  =  -0.75  X. 
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Figure  34.  Theoretical  and  measured  amplitude  of  Ez  at  three  distances  from 
the  array  for  x  =  0.75  X. 
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Figure  35.  Theoretical  and  measured  phase  of  Ez  at  three  distances  from  the 
array  for  x  =  0.75  X. 
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Figure  36.  Theoretical  and  measured  amplitude  of  Ez  at  three  distances  from 
the  array.  Mutual  coupling  effects  are  included  in  the  feed. 
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Figure  37.  Theoretical  and  measured  phase  of  Ez  at  three  distances  from  the 
array.  Mutual  coupling  effects  are  included  in  the  feed. 
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Figure  38.  Theoretical  and  measured  amplitude  of  Ez  at  the  boundaries 

(y  =  2.0  X,  3.5  X)  and  the  center  (y  =  2.75X)  of  the  more  distant 
test  volume. 
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Figure  39.  Theoretical  and  measured  phase  of  Ez  at  the  boundaries  (y  =  2.0  X, 
3.5  X)  and  the  center  of  the  more  distant  test  volume. 
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Figure  40.  Theoretical  and  measured  amplitude  of  Ez  at  the  boundaries  and  the 
center  of  the  more  distant  test  volume.  Mutual  coupling  effects 
are  included  in  the  feed. 
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Figure  41.  Theoretical  and  measured  phase  of  Ez  at  the  boundaries  and  the 
center  of  the  more  distant  test  volume.  Mutual  coupling  effects 
are  included  in  the  feed. 
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Figure  42.  Measured  amplitude  of  Ez  at  X/4  spacings  throughout  the  test 
volume. 
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Figure  43.  Measured  phase  of  E7  at  x/4  spacings  throughout  the  test  volume 
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Figure  44.  Measured  amplitude  of  Ez  at  A/4  spacings  throughout  the  test 
volume  for  x  =  0.75  A. 
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Figure  45.     Measured  phase  of  Ez  at   X/4  spacings  throughout  the  test   volume 
for  x  =  0.75  X. 
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Figure  46.  Measured  amplitude  of  Ez  at  X/4  spacings  throughout  the  test 
volume  for  x  =  -0.75  X. 
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Figure  47.     Measured  phase  of  Ez  at   X/4  spacings  throughout  the  test  volume 
for  x  =  -0.75  X. 
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Figure  48.  Measured  amplitude  of  Ez  at  x/4  spacings  throughout  the  test 
volume.  Mutual  coupling  effects  are  included  in  the  feed. 
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Figure  49.     Measured  phase  of  Ez  at   X/4  spacings   throughout  the  test   volume, 
Mutual    coupling  effects  are  included  in  the  feed. 
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Figure  50.  Measured  amplitude  of  Ez  at  A/4  spacings  throughout  the  more 
distant  test  volume  (2.6  A  <  y  <  3.5  A). 
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Figure  51.  Measured  phase  of  Ez  at  A/4  spacings  throughout  the  more  distant 
test  volume. 
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Figure  52.  Measured  amplitude  of  Ez  at  A/4  spacings  throughout  the  test 
volume.  Mutual  coupling  effects  are  included  in  the  feed. 
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Figure  53.  Measured  phase  of  Ez  at  A/4  spacings  throughout  the  test  volume, 
Mutual  coupling  effects  are  included  in  the  feed. 
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and  development  in  those  disciplines  of  the  physical  and  engineering  sciences  in  which  the  Bureau  is  active. 
These  include  physics,  chemistry,  engineering,  mathematics,  and  computer  sciences.  Papers  cover  a  broad 
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Handbooks — Recommended  codes  of  engineering  and  industrial  practice  (including  safety  codes)  developed  in 
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of  materials,  compiled  from  the  world's  literature  and  critically  evaluated.  Developed  under  a  worldwide  pro- 
gram coordinated  by  NBS  under  the  authority  of  the  National  Standard  Data  Act  (Public  Law  90-396). 
NOTE:  The  Journal  of  Physical  and  Chemical  Reference  Data  (JPCRD)  is  published  quarterly  for  NBS  by 
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Building  Science  Series — Disseminates  technical  information  developed  at  the  Bureau  on  building  materials, 
components,  systems,  and  whole  structures.  The  series  presents  research  results,  test  methods,  and  perfor- 
mance criteria  related  to  the  structural  and  environmental  functions  and  the  durability  and  safety 
characteristics  of  building  elements  and  systems. 

Technical  Notes — Studies  or  reports  which  are  complete  in  themselves  but  restrictive  in  their  treatment  of  a 
subject.  Analogous  to  monographs  but  not  so  comprehensive  in  scope  or  definitive  in  treatment  of  the  subject 
area.  Often  serve  as  a  vehicle  for  final  reports  of  work  performed  at  NBS  under  the  sponsorship  of  other 
government  agencies. 

Voluntary  Product  Standards — Developed  under  procedures  published  by  the  Department  of  Commerce  in 
Part  10,  Title  15,  of  the  Code  of  Federal  Regulations.  The  standards  establish  nationally  recognized  re- 
quirements for  products,  and  provide  all  concerned  interests  with  a  basis  for  common  understanding  of  the 
characteristics  of  the  products.  NBS  administers  this  program  as  a  supplement  to  the  activities  of  the  private 
sector  standardizing  organizations. 

Consumer  Information  Series — Practical  information,  based  on  NBS  research  and  experience,  covering  areas 
of  interest  to  the  consumer.  Easily  understandable  language  and  illustrations  provide  useful  background 
knowledge  for  shopping  in  today's  technological  marketplace. 
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